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ABSTRACT 

We derive a new age for the 7^ Velorum binary by comparing recent observations to 
our set of binary models. We find that it is very unlikely the stars have not interacted, 
which implies that previous estimates of the age from single-star models of 3.5 ± 0.4 
Myrs are incorrect. We prefer an older age of 5.5 ± 1 Myrs that agrees with the age 
of other lower mass stars in the Vela OB association. We also find that our favoured 
binary model shows the components of the binary have interacted in a Case B, post 
main- sequence, mass-transfer event. During the mass-transfer event, the envelopes 
of both components where radiative and therefore the damping of tidal forces are 
relatively weak. This explains why the binary is still eccentric after mass-transfer. 
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1 INTRODUCTION 

The binary 7^ Velorum contains the closest example of a 
Wolf-Rayet (WR) star to the Sun. It is a double-lined spec- 
troscopic binary and an extremely useful tool in understand- 
ing the evolution of Wolf-Rayet stars as the masses, lumi- 
nosities and m ass- loss rate can al l be determined. Recent 
observations bv lNorth et al] l|2007l ) have refined the param- 
eters of the binary orbit and the stellar components. A key 
parameter that was estimated was the lower limit age of 
bot h stars to be 3 . 5 ± .4 Myrs. 

IJeffries et al] (|2009l ) have took this lower limit of 3.5 
Myrs as the age and suggested that the system was one 
of the last to form in the Vela OB association, based on 
the observation that the other stars are estimated to be 
older than 5 Myrs. The age of 7^ Velorum was derived by 
fitting non-rotating single-star isochrones to the secondary. 
This method assumes that the stars have never interacted 
and the stars are not rotating. For an average rotation rate 
the age would inc r ease by about 25 percent to 4.4 Myrs 
iMevnet fc Maederl (I2OO0I '). Then assuming no binary inter- 
action initially seems reasonable as the system has an or- 
bital eccentricity of 0.3 and naively this means that the sys- 
tem never came into contact. Further consideration however 
finds that the orbital separation is small enough, ~ IOO-R0, 
that it is almost impossible for the stars not to interact 
through either Roche-Lobe Overflow (RLOF) or the more 
extreme Common-Envelope Evolution (CEE). We therefore 
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suggest that the two stars did interact in a timescale short 
enough that tidal circularization either did not occur or was 
not strong en ough to ach i eve a completely circular orbit. 
We note that iNorth et al.l (|2007l ') do state that such inter- 
action is possible and implies that their derived age should 
only be considered a lower limit. This was not considered by 
IJeffries et all (|2009l ) when they used their age in their study. 

We note that the observed parameters of 7^^ Velorum 
depend upon the assumed distance for the sys tem. In this 
work we have used the distance derived by iNorth et al.l 
(120071 ') from fitting the binary system of 336^7 parsecs. This 
is greater than the Hipparcos distance of 258^3^ parsecs 
and similar to the distance estimate from a single measure- 
ment by the V LTI that gave a dis tance to 7^^ Velorum of 
368^^3 parsecs l|Millour et al.ll2007l '). We have accounted for 
the uncertainty in the distance to the binary by including 
an error of ±0.1 dex in the luminosities of the binary com- 
ponents. 

There is a large amount of literature showing 
how binary interactions can change t he evolution of 
both s tars in a binar y llPod sia dlowski. Joss fc Hsul 
199^: IWellstein fc Langerl Il999l: IVanbeverenI l200ir 
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different sources of single-star e volution models an d they 
mostly agree with each other (e.g. lSmartt el al]l2009l '). there 
are no widely available binary equivalents. Construction of 
a set of binary models is not a simple process because rather 
than only needing to vary initial mass to obtain different 
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evolution, the initial binary separation and initial mass 
of the secondary must also be varied. This increases the 
computer time required to cover the full range of possibly 
evolutionary paths. 

Here we use the set of detailed b i nary evolution models 
created by lEldridee. Izzard fc ToutI ()2008h to determine a 
new age for •y'^ Velorum. First, considering the possible ages 
of the 9M0 WR star. Second, searching through our binary 
models for systems that agree with the masses and orbital 
separations of the system today. Third, the age of the sec- 
ondary star is considered. Afterwards, our results and the 
implications of the new age we derive are discussed. 

2 ESTIMATING THE AGE FROM STELLAR 
EVOLUTION MODELS 

The stellar models we used were calculated with 
the Cambridge STAR S cod e and were described in 
lEldridge. Izzard fc Tout! (|2008l '). Here, we restrict ourselves 
to use models with a metallicity mass fraction ot Z = 0.020, 
which we take to be Solar. These model s are all based on cir- 
cular orbits. However as discussed bv iHurlev. Tout fc Pols! 
(|2002l ) circular and eccentric orbits with the same semi-latus 
rectum should produce equivalent evolution. 

2.1 The age of the Wolf-Rayet star 

It is difHcult to accurately determine the age of a WR star. 
The star has lost at least half of its initial mass which could 
have been 20 to IOOMq. However it is possible to estimate 
an age range for a WR star from its current mass and WR 
su btype. The star is de signated as WRll in the catalogue 
of Ivan der Huchtl (|200ll ) and is identified as a WC star. This 
means it is a highly stripped star that has lost all its hydro- 
gen and most of its helium with the atmosphere dominated 
by carbon. We therefore search through our stellar mod- 
els and record the age of WC stars with masses between 8.4 
and 9.6M0 , the inferred mass of the WR star in 7^ Velorum. 
Here we consider the WR star alone and not the compan- 
ion or binary orbit. We assume a model is a WC star when 
there is no hydrogen present and {Xc/3 + Xq /4:)/Y ^ 0.03, 
where Xc, Xq and Y are the mass fra ction abundance of 
carbo n, oxygen and helium respectively (|Maeder fc Mevnetl 
1 19941 ). 

We find that our models indicate a range of possible 
ages from 3.3 to 6.4 Myrs for WC stars with masses simi- 
lar to that of the WR star in j'^ Velorum. Ou r lower esti- 
mate of 3.4 Myrs agrees with the age derived bv lNorth et al.l 
l|2007f ) for the O star, although this was with older single star 
isochrones. The youngest WC stars are from stars with ini- 
tial masses above 6OM0 and the older WC stars are from 
stars with an initial mass of 3OM0. star. Lower initial masses 
are preferred by the initial mass function thus initially less 
massive, and therefore older, stars are the preferred when 
we estimate the age of WRll. 

2.2 The possible initial parameters for the 7^ 
Velorum binary 

We have searched through our binary models to find those 
with a reasonable match for the 7^ Velorum system. For a 



certain binary model to match, we require that the following 
are all true at some point during the lifetime of the binary. 

• The WR star is of WC type. 

• The mass of the WR star is between 8.4 and 9.6M0. 

• The luminosity of the WR star is between log(L/L0) = 
5.1 and 5.3. 

• The secondary mass is between 27 and 3OM0. 

• The orbital separation is between log(a/J?0) = 1.8 and 
2.3. 

• Any interaction occurs after the main sequence (Case 
B mass-transfer). 

As already discussed the binary is eccentric and all our 
models have circular orbits. We therefore assume that, for 
our binary systems to fit, it must have a radius somewhere in 
the range of separations of the observed binary. The effect 
of eccentricity woul d be for mass-transfer to begin earlier 
(|Church et al.l [20091 ) and thus the initial separation would 
be greater for more eccentric systems to obtain the same 
evolution as we discuss here. 

We choose to only consider systems with mass-transfer 
after the main sequence (Case B). This is because when 
mass-transfer events occur on the main sequence (Case A) 
they happen on a nuclear timescale and the tidal forces 
would have more time to circularise the orbit. Post main- 
sequence mass-transfer occurs on a thermal timescale and 
therefore tidal forces have less time to circularise the orbit 
and the resulting binary remains eccentric. 

We do not attempt to model the secondary here be- 
yond matching its mass. The secondary does lose mass by 
stellar winds and mass is transferred from the primary to 
the secondary. The accretion rate onto the secondary is lim- 
ited to A/2 / Tthcrmai where Tthcrmai IS the sccoudaries thermal 
timescale. If the mass-transfer rate is greater than this value 
the excess mass is lost from the system. We do not include 
angular momentum transfer and thermohaline mixing in 
this model. This is because they produce complex affects on 
the secondary, especially angular momentum transfer which 
will alter the amount of rotationally induced mixing. This 
makes the evolutionary outcome of the secondary uncertain 
ICantiello et al.ll2007l : [Stancliffe fc Eldridgell2009l ). We con- 
sider the secondary in more detail in Section 2.3. 

We list our models that match 7^^ Velorum in Table [1] 
The range of initial primary masses and initial separations 
are narrow. The initial separations are similar to the ob- 
served separation. This indicates that the widening of the 
system due to stellar-wind mass-loss is countered by the 
mass-transfer event tightening the binary. 

For each of the systems a binary interaction must have 
occurred. Therefore we might simply assume that 7^ Velo- 
rum should be a circular orbit rather than eccentric. It is 
worth asking, how unusu al it is for a Woff-Ray et binary to 
be eccentric. Figure 3 in Ivan der Huchtl l|200ll ) shows that 
in general binaries with periods of 30 days or less tend to be 
circular. However there are a number of systems with peri- 
ods between 30 and 100 days that have eccentricities in the 
range of 0.3 to 0.6. This is because binaries with periods be- 
low 30 days experience Case A mass-transfer and therefore 
their tidal forces have a long time to circularise the orbit. 
The binaries with longer periods experience Case B mass- 
transfer, that proceeds at shorter thermal timescales giving 
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Table 1. The systems in our grid of binary models that match 7- Velorum. We list their initial parameters and the range of the parameters 
when they match the observed details of 7- Velorum. 
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log(a/ij0) 


Mi/Mq 


M2/MQ 


/Myrs 


log(a/ijQ) 


[C/Hc 


■]WR 


[O/Hc 


■]wR 




35 


31.5 


2 


9.6 - 8.4 


30.0 - 29.8 


5.51 - 5.59 


1.81 - 1.82 


0.03 - 


0.25 


0.001 - 


0.031 


5.24 - 5.16 


35 


31.5 


2.25 


9.6 - 8.4 


29.8 - 29.6 


5.51 - 5.59 


2.06 - 2.08 


0.05 - 


0.30 


0.003 - 


0.045 


5.25 - 5.17 


40 


20 


2 


9.6 - 8.4 


30.0 - 30.0 


5.01 - 5.13 


1.89 - 1.91 


0.27 - 


0.55 


0.036 - 


0.153 


5.26 - 5.20 


40 


28 


2.25 


9.6 - 8.4 


27.5 - 27.5 


5.01 - 5.14 


2.07 - 2.09 


0.28 - 


0.57 


0.040 - 


0.168 


5.26 - 5.20 



less time for tidal forces to circularise the orbit. These facts 
suggests that ■y'^ Velorum is not unusual. 

From the above, the assumption that binary mass- 
transfer events must always circularise the orbit is not true. 
We have estimated a circularization timescale for the mod- 
els in Table [T] and compared this to the time for the CEE. 
The key to determine the timescale of circularization in a bi- 
nary is the damping mechanism for the tides created in the 
star. In general stars with a convective envelope have shorter 
timescale than those with radiative envelopes by around tw o 
orders of magnitude l|Zahnlll977l : lHurlev. Tout fc Polsll2002l ). 
We find that in our 35Mq model the star is mostly radiative 
with a radius of 90Rq at the time of interaction. The model 
does have a small convective shell just outside the helium 
core between 1 to 4i?0 but it is too deep within the star to 
effectively damp tides. 

Using the details giv en in iHurlev. Tout fc Pols! (I2OO2I ') 
and Ide Mink et all ||2009D . we have estimated the circular- 
ization timescale of the binary to be a few thousand years, 
similar to the time of the interaction in our model. The 
two timescales are similar in order of magnitude, but be- 
cause both are uncertain we conclude it is unlikely that 
complete orbit circularization is possible. In comparison we 
find that for initially tighter binaries mass-transfer occurs on 
the main-sequence via RLOF only over a timescale of > 0.1 
Myrs. Therefore there would be plenty of time for complete 
circularization to occur. 

One final detail to consider in Table [T] are the ratios 
of surface carbon, oxygen and h elium a bundances by num- 
ber from our models. iDc Marc ho et al.l 1^)00) derived the 
abundance ratios for 7^ Velorum from different model atmo- 
sphere codes. They found that [C/He] took values of 0.15 to 
0.06 and [O/He] took values of 0.03 to 0.01 respectively. The 
binaries in Table [T] agree with these inferred compositions. 
However looking at the models in detail when we match the 
[C/He] ratio we underestimate [O/He] by a factor of 3 indi- 
cating the rotation may play some role in the evolution of 
the primary star during the WR phase. 

2.3 Modelling the secondary 

Modelling the secondary is not an easy task. If the two 
stars have interacted in some way it is difficult to know 
exactly what effect this would have had on the secondary 
l|Cantiello et al.ll2007l : IStancliffe fc Eldridgd [20091 ) . Here we 
restrict ourselves to study the effect of rotation on the sec- 
ondary alone. We see from Table[l]that there is a wide range 
of initial masses possible for the secondary. Some possible 
progenitors accrete a large amount of mass to reach that ob- 
served today. These stars with lower initial mass, M2 — 20 



and 28M0, are less evolved at the time of mass-transfer 
and thus have lower luminosity and/or higher surface tem- 
perature at the time when the primary star matches the 
observations. 

Here we concentrate on the secondary star from our 
favoured model with an initial mass of 31.5M0 which ac- 
cretes very little mass. We find that this star has the same 
mass as observed for the secondary at the same time as the 
primary is a WC star with a mass of 9M0 and fulfills our 
criterion outlined in Section 2.2. However as shown in Figure 
[T]a basic model only just agrees with the observed luminos- 
ity and radius and 7^ Velorum would need to be closer than 
300 parsecs. This disag rees with the distanc es measured by 
iNorth et all (|2007l ') and lMillour et al.l l|2007h . If we consider 
that the star can only be more distant th an we have as- 
sumed as indicated bv lMillour et al. I 12007') and thus more 
luminous it suggests that our basic model cannot reproduce 
the observed secondary radius. In Figure [5] we see that the 
basic model agrees with the luminosity and mass simultane- 
ously but only at the very end of its main-sequence lifetime. 
Using Figure [3] we estimate an age from the basic secondary 
model of 5.6^?;3 Myrs. 

Rotation is a complicated physical process in stars and 
our code does not currently include it. The general effect 
of rotation is to induce mixing in a star that increases the 
time spent on the main-sequence by mixing nuclear pro- 
cessed material to the surface of the star and new hydrogen 
into the core. To simulate this effect we have included an ar- 
bitrary small amount of extra-mixing through the radiative 
zone of our model. We see in Figure [T] that with this mix- 
ing the star achieves a higher luminosity at a given radius 
and eventually agrees with the observed luminosity and ra- 
dius of th e secondary. Comparing to models that do include 
rotation (|Mevnet fc Maeden l2005h the evolution is similar 
to that provided by reasonable rotation velocities around 
200km s"^ 

We see in Figure [2] that the model with enhanced mix- 
ing also agrees with the observed mass and luminosity of the 
secondary and is in the middle of its main-sequence evolu- 
tion. Figure [3] allows us to estimate the age of the secondary 
star to 5.4^j'2 Myrs. This is consistent with the range of 
ages we determined for the age of the WC star and binary 
system above. 



3 DISCUSSION AND CONCLUSIONS 

We have used three different methods to estimate the ages 
of 7^- Velorum. These were estimating the age of a 9M0 
WC star, modelling the binary orbit and components simul- 
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Figure 1. The luminosity-radius relationship for two different 
stellar models with initial masses of 31.5M0. The black line indi- 
cates our normal stellar model while the red dash-dotted line is for 
a stellar model where a small amount of mixing is allowed within 
the radiative zones, simulating the effect of rotation on the evo- 
lution. The dashed line indicates the parameters of the secondary 
in the 7-Velorum binary with the grey shaded region indicating 
the error in the luminosity due to the uncertain distance. 



Figure 3. Similar to Figure [T] but the relationship between 
and luminosity of our stellar models.. 
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Figure 2. Similar to Figure[T]but the relationship between stellar 
mass and luminosity of our stellar models. 



taneously and considering the age of the best fitting initial 
secondary mass with and without enhanced mixing. We find 
that the age must be older than 3.5 Myrs and by combin- 
ing different ages we estim ate that it is 5.5 ± 1 Myrs. The 
2 Myrs difference is because iNorth et al.l |20o3) used single- 
star non-rotating isochrones to determine a lower limit for 
the secondary star age. This work demonstrates the neces- 
sity for publicly available grids of binary star models to be 
created. The use of single-star models can give misleading 
results that lead to incorrect conclusions. 

Our age estimate is in better agreement with 
ages derived for the surrounding lower-mass stars from 



Ijeffries et al. 

I 120091). It suggests that all the stars in the 
Vela OB2 association formed at a similar time. 

The older age estimate may also help explain the non- 
detection of 1.8 MeV photon emission from 7 Velorum. This 
emission is due to the radioactive decay of '^^ Al formed dur- 
ing core hydrogen burning by proton capture on '^^Mg. ^^Al 
ha s a half-life of 0.75 Myrs and 7 Velorum is predicted 
by iMowlavi fc Mevnetl l|2006l ) to have detectable emission. 
Currently no such emission has been detected. We suggest 
there are two factors reducing the amount of Al below the 
detectable abundance. First in our models of ■y'^ Velorum 
we find the time difference between the end of core hydro- 
gen burning and our age estimate is 0.85 Myrs thus less 
than half of '^^ Al created remains. Second the lower initial 
mass of the WR star may mean less ^^Al is created dur- 
ing the main-sequence as predicted for the 60Mq star by 
iMowlavi fc Mevnetl (|2006l ). 

The fact that the binary is still eccentric does not rule 
out a binary mass-transfer event. Our binary models indi- 
cate that some post main-sequence (Case B) mass-transfer 
events occur so rapidly that the tides formed in the mostly 
radiative envelope are unable to circularise the orbit. On 
the other hand, WR binaries with periods below 30 days 
experiences Case A mass-transfer on the main-sequence and 
the tides do have time to circularise the orbit. This places 
a limit on the circularization timescale of tides in radiative 
envelopes to a few thousand years, the time required for the 
hydrogen envelope to be removed in a mass-transfer event. 

Our major remaining uncertainty is the evolution of the 
secondary. We find that it is possible to explain its current 
evolutionary state if initially it was a 31.5Af0 star with a 
rotation rate of ~ 200km s~^. Less massive stars that ac- 
crete a large amount of mass never match the radius and 
luminosity of the secondary when the primary star is a WC 
star of 9M0 . Future study of this system may provide more 
clues to the effect of rotational mixing and of mass-transfer 
on the secondary stars in binary systems. 
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